Phase coherence and fragmentation in weakly interacting bosonic gases 
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We present a theory of measurement-induced interference for weakly interacting Bose-Einstein 
condensed (BEC) gases. The many-body state resulting from the evolution of an initial fragmented 
(Fock) state can be approximated as a continuous superposition of Gross-Pitaevskii (GP) states; 
the measurement breaks the initial phase symmetry, producing a distribution pattern corresponding 
to only one of the GP solutions. We discuss also analytically solvable models, such as two-mode 
on-chip adiabatic recombination and soliton generation in quasi one-dimensional condensates. 
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A long-standing fundamental problem in theoretical 
physics is to understand how relative phases are estab- 
lished between superfluids that have never been in con- 
tact with each other In the case of atomic Bose- 
Einstein condensates, the first experiment of the sort 
was done already a decade ago [|J ; strangely though, the 
data could be reproduced simply by assuming a phase 
relation between the initial condensates and using the 
time-dependent GP equation It soon became ac- 
knowledged, following a number of elegant proofs for the 
noninteracting case [J, that the U(l) phase symme- 
try is broken by the measurement process itself. Since at 
the time of the measurement the density of the expand- 
ing condensates is relatively low, it is assumed that the 
situation is for all practicalpurposes equivalent to the 
noninteracting case. Since 2], better-controlled experi- 
ments in various setups (optical lattices [6( , on-chip split 
condensates [7j, nondestructive measurements |8||) have 
supported this view. 

However, in the interacting case the situation has 
proven to be more complicated: on one hand, it is known 
that interaction within each cloud tends to delocalize the 
phase between two condensates, leading to squeezing and 
to other effects similar to the Josephson effect in super- 
conductors @. But BEC interference experiments re- 
quire overlap of the atomic clouds, therefore the physics 
at work might well be different. If atoms originating from 
the initially independent clouds interact, it is not even 
clear in what sense, at the measurement time, we can talk 
about condensates that have not seen each other. Since 
fragmented states are not robust to certain classes of per- 
turbations, it could well happen that the interaction it- 
self would lead to phase localization [l(|. Moreover, for 
repulsive interactions a typically phase-coherent ground 
state is energetically favored H,Hl| ■ Very recently, some 
authors have calculated, using standard many-body tech- 
niques pjl the average of the density operator on evolv- 
ing fragmented states and found out that in the interact- 
ing case the average of the density operator could display 
ripples which look somewhat similar to the fringes ob- 



served in atomic interference experiments. It is yet not 
clear if a phase-coherent condensate is formed due to in- 
teractions before the measurement starts [T2j. Could it 
be then the case that what is detected in BEC interfer- 
ence experiments with interaction is in fact the density 
and not higher-order correlations? 

In this paper, we present a theoretical approach that 
takes into account both the effects of interaction and that 
of measurement. We show that, also for weakly interact- 
ing bosons, the broken-symmetry approach to this prob- 
lem H isiustifiable, and there is no contradiction with 
results [12j showing ripples in the average particle den- 
sity. 

The Hamiltonian of a system of bosons with repulsive 
two-body interaction is 



H 



dr^>' (r)iJoV>( r ) 



drtjj' (r)^^(r)^(r)^j(r), 



(1) 

with g = A-KWa/m > 0, and H = {-h /2ra)V 2 + V cxt (r) 
the Hamiltonian of a single particle moving in an exter- 
nal potential V C xt- We start by presenting the general 
argument for evolution and measurement; to gain further 
insight, we examine later in the paper particular cases in 
which the GP equation can be solved analytically. 

Evolution We first show that, starting from an initial 
fragmented state | gnd L ) | gnd fl ), the time-dependent 
many-body state can be approximated as 



exp[--Ht] | gnd L ) | gndf 
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2^ 



N, 



(2) 



'Electronic address: paraoanu@cc.hut.fi 



where c = 2 w / 2 (iV/2)!/v / M @, | %(t)) N = 

, r i N 

(VM)" 1 /d 3 r$ v (r,f)^( r ) | o), and ^(r,*) will 
be identified as the result of evolving the initial state 

$ v (r,0) = (1/V2) (<Mr)e-^/ 2 + <Mr)e^/ 2 ) by the 
time-dependent GP equation. We will proceed in a sys- 
tematic way, using a perturbative expansion technique 
developed in a different context [H, [3]. We first con- 
struct, at any time t and for a given phase ip, a projection 
operator R v (t) outside the subspace spanned by the ket 
| $ v (t)), R v (t) = 1- | |, and we split the 



2 



field operator tp(r') = $ v (r', t) + \ v (r', £), where 

X v (r',f) - / dTR v {v\v-t)^{v). (3) 



Next, we linearize the Hamiltonian fl} around $ v and, 
after relatively lengthy but straightforward calculations 
we find that in the Heisenberg picture (denoted by 
a superscript (W)), the evolution of the phonon operator 



X^ h \y' ,t) takes the form 

*ftj^#V.*) = / drR v (v,v;t) 



-in—%(r;t) 

(H) 



+H %(v,t)+g$Ur,t)$Jr,t)%(r,t) , (4) 



and therefore 



1 ^ { ^' 1 $v>(<) ^'= v rV2 / H(r'M)[-4^( r ' f ) +iTo4 v (r,t) + . 9 A | %(r,t) | 2 * v (r,t)] (5) 
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The last equality is written up to an (irrelevant) phase 
factor. To obtain it, we notice that, as an immediate 
application of the Cauchy-Bunyakovski-Schwarz inequal- 
ity in L 2 Banach spaces, | ($ v '(t) \ ® v {t)) | reaches 
a maximum value of 1 for iff = ip. The quantity 
1 7v (®cp'(~t) I ®<p(t))N I will then be very close to zero 
for ip' ^ ip, and strongly peaked to 1 if ip' = ip, so we can 
write 



I (V I #„)(*) \ v >= v (v?-<p) 



, N d 2 

exp -y I ^72 I (V I \<p>=v\ if' - V) 2 



The selfconsistency condition for the validity of the 
expansion Eq. ([2]) is that at any time the states $ v (r, t) 
remain macroscopically occupied, 



^^(VWI^lr'OI^W^ (6) 
» / Vw-i^'W I X v (r',i) | M*)>w (7) 



or in other words the quantum depletion of a state tp 
should not aquire contributions of the order of the corre- 
sponding order parameter during evolution. According to 
Eq. |5]), this is realized if $ v (r, t) is a solution of the time- 
dependent GP equation, ih(d/dt)$ v (r,t) — Ho$> v (r,t) + 
gN | $ v (r,t) | 2 $ v (r,t) with initial condition $ v (r,0) = 
^ (^(r^"^/ 2 + (j) R {v)e^/ 2 ). This is because the evo- 
lution Hamiltonian has vanishigly small matrix elements 
between states | $ v (t))N with different phases ip; there- 
fore each of these states is evolved independently. The 



limits of validity of this formalism are identical to those 
of the time-dependent Gross-Pitaevskii equation for a 
macroscopic occupation of <&<p, namely that dynamical 
instabilities do not develop (the Bogoliubov-de Gennes 
equations should not have imaginary eigenvalues). We 
also note that, using the orthonormality property de- 
scribed above, the same result can be obtained by min- 
imizing the action JJ" dt[(l/2)(ih(d/dt)) + c.c. — (H)] on 
the state Eq.©. 

Measurement Suppose now that we have a series of 
n ^> 1 detections leading to coordinates ri, r 2 , r„; the 
resulting wavefunction is then 



/n 
dip 

3=1 



(8) 



Given a certain sequence ri, T2, r„, there are two pos- 
sibilities: either there is a wavefunction $ v (r, t) whose 
shape it approximates best - and in this case this wave- 
function is selected - or there is no such wavefunction and 
in this case the probability of this particular sequence is 
very small. To prove this statement, we proceed by math- 
ematical induction. Let us divide the detection space into 
small bins of volume v each, centered at discrete positions 
rj. A particle can be detected in either one of the bins i. 
Suppose that after n detections we have rii atoms at r,, 
resulting in the wavefunction 



^n^^'*)riM*)> 



N—n • 



(9) 



We now show that the maximum absolute value of the 
coefficient Yli [®ip{ r h t)) nz corresponds to a tp for which 
rii/(nv) =| <&(p(Yi,t) | 2 , in other words the "histogram" 
of rii 's is mimicked by the Born rule for the wavefunction 
corresponding to (p. We use the method of Lagrange 
multipliers, with the normalization | $ v (rj,t) | 2 v = 



3 



1 as a constraint; we then have 



d_ 

dip 



A solution of this equation is obtained for | $^(rj,t) | 2 
/ | $$(rj,t) \ 2 = rii/rij; supplemented with the nor- 
malization condition, this produces the announced result 
rii = nv | $$(ri,t) | 2 . 

This means that the next measurement event will tend 
to reduce even more the probability amplitude of states 
which differ significantly from | $^(r i; t) |, since these 
states have anyway a smaller probability amplitude to 
start with. 

We will now show that, as n increases, the probabil- 
ity amplitudes ]X | <& y (rj,i) \ Ui peak strongly around a 
value corresponding to (p. We make a Taylor expansion 
around the maximum value, 



where [Y] is a column vector with transpose \Y] t — 
(..., $m(rj, t) — $p(rj, i), ...), and VW is the Hessian ma- 
trix with elements M.a = (rii — l)?w]X \ \ ni 
and Mij = ^/rhnjH, j *tp(r l; i) |™ ! (i 7^ j). Using the 
constraint imposed by normalization, we find 
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Finally, averaging over many shots eliminates the 
symmetry-breaking effect, but in general it does not flat- 
ten out all the ripples since those related to the creation 
of collective excitations (an effect due only to interac- 
tions) will survive (see also the model below); thus there 
is no contradiction with the findings of 

On-chip condensates As a toy model in which the cal- 
culations can be done analytically, we consider the case 
of two condensates launched in the lowest modes 4>i and 
4> r of two on-chip atomic waveguides (see e.g. Fig. f in 
[15l|). In this case, it is possible to make an (adiabatic) 
two-level approximation and the problem can be solved 
analytically. In the conditions specified in only the 
modes <j>\ — (4>i+4> r )/2 and (f>2 — {<j>i—<j> r )/2 (correspond- 
ing to two detection channels similar to Q) participate 
in the dynamics. 

Consider now what happens after the system (ini- 
tially in a fragmented state) has evolved for some time; 
the complex amplitude probabilities corresponding to 
each mode are of modulus ci = y(T— y)/2, c-i = 



\/ {1 + y)/2 where y is a parameter that depends on time, 
on the initial phase difference ip, and on the nonlinearity 
(see [l5l | for explicit notations). Using the notations of 
[15| . we now imagine, in the spirit of [5j, that we make k\ 
detections in the 4>i channel and &2 detections in the (f>2 
channel: the resulting (unnormalized) many-body func- 
tion becomes 



2 71 



(10) 



The maximum value of the quantity c^c^ 2 is reached 
for y — (&2 — ki)/k, where k — k\ + k2] expand- 
ing around the maximum value, we get c 1 1 c 2 2 ~ 

(h/kjf^ 2 (k 2 /k) k2/2 exp [-(fc 3 /i6M 2 )(y - yf] ■ 

After a large enough number of measurements, k 3> 1, 
the value of y tends to localize around y. Indeed, the 
worst-case scenario for the localization of y corresponds 
to a minimum value for the factor k 3 /kik2, which (with 
the restriction k\ + fc 2 = k), is 4fc, therefore increasing 
with the number of measurements. We then find 



ki k 
CV C 



fcl+1 



-f) 6(3, -y). (11) 



Now, if one wants to average over many single-shot 
measurements, the average number of particles detected 
say in channel 1 will be N ^ (dp /2-k)c\ 7^ N/2 which is 
the two-mode analogue of a density ripple due to inter- 
actions (e.g. for asymptotically large times f -> 00, this 
average is N, while the corresponding average for c\ is 
zero; in the noniteracting limit both averages would be 
N/2). 

Also, in the particular case of zero interaction, we 
recover the same results as in 5], using the substitu- 
tion y = cos tp, namely localization of the phase near 
<p G {arccosj/, 2tt — arccosy}, 



c l c 2 



7T / k\ 



6((p-<p). 



(12) 



Quasi one- dimensional condensates In quasi one- 
dimensional infinite uniform gases, the GP equation has 
two solitonic solutions propagating in the +z and — z di- 
rections respectively with speed v s — ^/2jicos(p/2. Nu- 
merical simulations (l6l.[l7j show that, when recombining 
two such quasi-condensates with a previously established 
phase difference between them, a train of solitons is in 
fact formed. Consequently, when two fragmented quasi- 
condensates interfere, our theory predicts the appearance 
of solitonic trains propagating (randomly) in the ± direc- 
tions. 

To get a better understanding of the left- right- local- 
ization process, let us simplify drastically the real phys- 
ical situation and consider that the initial density of the 
two components is such that, upon connection, only two 
solitons with notches at z — y / 2/iicos^/2 > and 
— zq = \f2~jlt cos(-7r — ip/2) < 0, corresponding to phase 



4 



differences of < ip < tt and 2n — tp, are preferentially 
formed (for other phases the propagating ripples will be 
too small to be measurable). Then the relevant part of 
the many-body wavefunction at t, when the detection 
is performed, has the generic form of a Schrodinger cat 
state, 



*} =1 $H 



N- 



N ■ 



(13) 



We show below that the symmetry ip — > 2n — ip of this 
state is broken by the measurement process. Absorption 
processes resulting in atoms detected far enough from 
the notches (where the wavefunctions $± differ only by 
phase factors) do not change significantly the probabilis- 
tic weight of the states | $+) and | $+), therefore they 
do not discriminate between the two solitons. Without 
loss of generality and to simplify the proof, we will dis- 
cuss only processes involving detection of atoms in either 
of the two notches at ±z . Consider a sequence of detec- 
tions, say {(+Zq), (— zq), (—Zq)....}, and define f s as the 
sign of the coordinate of the s atom detected. The prob- 
ability V {f} of a sequence {/} = {f 1} f 2 , f 3 , f n } can 
be calculated by applying N times the operator ip(z); we 
find 



V 



{/} 



n 

fe=2 



inti^-/.(/^o) i 2 

2E a =±intl 1$ a(/^0) 



(14) 



The sum of probabilities over all possible 2 N sequences 
is Ylf^f = !• Eq. ([H]) has the mathematical prop- 
erty that (for large N) Vf flattens to zero for sequences 
which contain comparable amounts of "+" and "-" de- 
tections, so it tends to favor one or the other of the states 
| <fr±). In other words, the final result of an experiment 
will be a distribution which will reflect the shape of ei- 
ther one of these states. This can be seen immediately 
for the particular case of approximately dark solitons, 
$ + (zo) = 0, $_(— z ) = 0, for which a single detec- 
tion event is enough to discriminate between the left- 



and right- propagating solitons. Indeed, after the first 
detection, the resulting (unnormalized) manybody wave- 
function is ${Z ) I *) = VN[<& + {z ) I $ + ) A r_ 1 +$_(z ) | 

$_)Ar_ 1 ] w| $_)tv-i if the atom is detected at zq, and 

$(-zo) I *) « VN[t> + (-z ) | $+)jv-i +$-(-*>) I 
3>_)jv_x] fa | $ + )at_ 1 if the atom is detected at — z . In 
the general case, our claim can be shown by mathemat- 
ical induction. Given a sequence / of n measurements, 
out of which n+ have been detections at +zq and rt_ at 
—zo, we can associate, according to Eq. (fT4|) . to the n+ 1 
detection event the probabilities 



P{/,±} = ^{/} 



;±2(n+-n_) 



(15) 



where f =| $_ (-z )/$_ (z ) \ = 
Then V{f,+} > if n+ > n 

n 



S_(-*o)/*-(*o) |< 1. 
and V {fM < V {f ^ } if 
< n_. This shows that the initial difference between 
n + and n_ will increase exponentially fast under subse- 
quent detections, which proves our point. As an example, 
an alternating sequence { + ,— ,+ ,— ,...,—} of n=even de- 
tections will have a probability 2~™/ 2 £™(l+£ 2 )~™/ 2 , while 
for a constant sequence { + , +, +, +} the correspond- 
ing probablility is + ^ k ~^]~ 1 ] one can imme- 
diately check that the ratio between the first probability 
and the second goes indeed very fast to zero as the num- 
ber of detections n increases. 

In conclusion, we have shown that the evolution of a 
fragmented-state BEC yields a many-body state which is 
a continuous superposition of GP states, embedding the 
effects of interaction and indistinguishability. An analy- 
sis of the measurement process reveals that this many- 
body state collapses onto only one of the GP states in 
the superposition. 
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